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An  Analysis  of,  and  Some 
Observations  on,  Dynamic 
Fracture  in  an  Impact  Test 
Specimen 

Numerical  simulations  of  crack-propagation  histories  in  four  cases  of  dynamic  tear 
test  experim  ents  on  4340  steel  are  performed.  The  influence  of  the  loss  of  contact  of 
the  specimen  at  various  times  with  either  the  supports  or  the  tup  or  both  is  critically 
examined.  In  each  case,  the  variation  of  the  dynamic  K-factor,  for  the  simulated 
crack-propagation  history,  is  directly  computed.  The  results  are  discussed  in  the 
light  of  current  controversies  surrounding  the  dynamic  fracture  toughness 
properties  governing  crack-propagation  under  impact  loading.  It  is  concluded  that 
these  controversies  may  not  be  fully  warranted. 


ABSTRACT 

-  Numerical  simulations  of  crack-propagation  his¬ 
tories  in  four  cases  of  dynamic  tear  test  experiments 
on  4340  steel  are  performed.  The  influence  of  the 
loss  of  contact  of  the  specimen  at  various  times  with 
either  the  supports  or  the  tup  or  both  is  critically 
examined.  In  each  case,  the  variation  of  the  dynamic 
K-factor,  for  the  simulated  crack-propagation  history , 
is  directly  computed.  The  results  are  discussed  lit 
the  light  of  current  controversies  surrounding  the 
dynamic  fracture  toughness  properties  governing  crack- 
propagation  under  Impact  loading.  It  Is  concluded 
that  these  controversies  may  not  be  fully  warranted. 

INTRODUCTION 

Until  recently,  for  stltuations  governed  by 
small-scale  yielding,  it  was  thought  that  the  govern¬ 
ing  criterion  for  elastodynamic  crack  propagation  un¬ 
der  Mode  I  plane  strain  conditions  can  be  written  as: 
Kj(v,t)  *  'rl*ere  Kjpfv)  *•*  the  velocity-depen¬ 

dent  fracture  toughness  of  the  material,  which  was 
thought  to  be  a  "reasonable"  geometry- Independent 
material  property.  This  hypothesis  appejred  to  have 
been  validated  In  several  studies  related  to  dynamic 
crack-propagation  Initiated  under  quasl-statfc  loading. 
In  the  analysis  of  such  cases,  botli  "generation"  and 
"propagation"  calculations  were  employed:  In  the 
former  calculation,  the  experimentally  measured  crack- 
propagation  history  was  simulated  to  find  the  stress- 
intensity  factor  or  the  velocity-dependent  fracture 
toughness:  the  latter  calculation  was  used  in  either 

of  the  two-wavs  (l)  based  on  a  given  K,_  versus  v 
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relation  to  find  the  crack  propagation  history,  or 
(11)  to  find  the  best  versus  v  relation,  the  cal¬ 

culated  crack-propagation  history  corresponding  to 
which,  agreed  best  with  the  experiment.  The  remarkable 
success  of  these  calculations  appeared  to  indicate 
that  the  prediction  of  dynamic  crack-propagation  and 
possible  arrest  under  general  loading  conditions  may 
be  well  within  the  grasp  of  current  art  of  computa¬ 
tional  mechanics. 

Recently  some  work  has  appeared,  however, 
that  seemed  to  cast  doubt  on  the  concept  of  dynamic 
fracture  toughness  that  Is  Independent  of  the  rate 
of  applied  loading.  In  Ref.(l),  Kanninen  et  al  re-  . 
ported  experimental  and  numerical  results  on  dyanamic 
tear  test  specimens  of  4340  steel,  a  high-strength,- 
rate-insensitive  material.  In  these  experiments, 
crack-propagation  was  initiated  from  notches,  of 
varying  degree  of  "bluntncss",  under  Impact  as  well 
as  quasl-scatlc  loading.  The  dynamic  fracture-tough¬ 
ness  was  attempted  to  be  inferred  from  the  energy 
measured  to  be  absorbed  In  the  tear  test.  A  series 
of  "propagation  type"  linear  elasto-dynamic  analyses, 
using  hypothetical  K  values,  were  performed.  A 
surprising  finding  of  [l|  was  that  the  dynamic  frac¬ 
ture  toughness  governing  crack-propagation  initiated 
from  a  blunted  crack-tip  under  impact  loading  may  bo 
significantly  higher  (roughly  170)  than  when  crack- 
growth  Is  Initiated  r|uasl-statlcallv  ( rough Iv  about 

65MNm'U5). 

Thu  primary  objective  ot  the  present  paper 
Is  an  attempt  to  analyze  the  data  presented  in  (U 
and  to  examine  the  results  in  the  light  ot  the  con¬ 
clusions  presented  in  [  1 1 .  In  addition  to  the  blunt- 
notch  specimen  data  reported  in  f  L  | ,  analysis  is  per¬ 
formed  also  ot  the  data  for  .1  fatigue  prc-c racked 
specimen  supplied  to  the  authors  bv  the  Buttelle 
Columbus  Laboratories  in  January  1981.  In  the  pres¬ 
ent  paper,  in  contrast  to  those  In  (1],  "’eneration* 
type  calculations  are  employed,  i.e.,  the  experimen¬ 
tal  data  tor  'rack-velocity  versus  time  historv  t* 
simulated  In  1  l Init  *  element  program  to  deterram- 


directly  the  stress-intensity  factor  variation  with 
time.  The  employed  finite  element  method  is  the 
"moving-singularity"  procedure  reported  earlier  by 
a  the  authors  [2,  3].  In  the  present  analysis,  care- 
'  ful  attention  is  paid  to  the  boundary  conditions  on 
the  specimen,  especially  the  loss  of  contact  of  the 
specimen  at  various  times  with  either  the  supports, 
or  the  tup,  or  both.  Four  different  cases  of  experi¬ 
mental  specimen  are  analyzed.  In  each  case,  the  In¬ 
put  energy,  kinetic  energy,  elastic  energy,  and  frac¬ 
ture  energy  variations  are  computed.  It  is  noted 
that  each  of  these  energy  quantities  is  computed  dt- 
rect Iv  in  the  present  procedure.  That  these  energy 
quantities  "blance  out"  is  nothing  but  an  a  posteriori 
check  on  the  present  calculations.  It  should  thus  be 
noted  that,  in  contrast  to  the  "propagation"  calcu¬ 
lations  in  [1],  fracture  energy  is  not  inferred  from 
an  energy  balance,  but  directly  computed  in  the  pres¬ 
ent  "generation”  calculation. 

Detailed  results  are  presented  for  each  of  the 
four  cases  analyzed.  These  results  are  analyzed  to 
arrive  at  some  "plausible"  conclusions  which  appear 
to  be  at  variance  with  the  conclusions  presented  in 
111. 

ANALYSIS 

The  test  specimen  geometry  is  indicated  in  Fig. 

1,  along  with  the  finite  element  mesh  employed  in  the 
modeled  portion  of  the  specimen.  Points  L  and  S  in 
Fig.  1  represent,  respectively,  the  loading  and  sup¬ 
port  poincs.  Sixty-cwo  8-noded  isoparametric  ele¬ 
ments  and  one  moving  singularity  element  are  used. 

The  specimen  geometry  indicated  in  Fig.  1  corresponds 
to  chat  reported  in  [1],  and  a  plane-strain  condition 
is  Invoked  in  the  present  two-dimensional  analysis. 

In  simulating  the  experiments  (1),  the  following  in¬ 
itial  conditions  are  used  in  the  present  analysis: 
at  time  t-0,  velocity  u^-6. 88m/ sec.  The  tup  displace¬ 
ments  are  calculated  by  u^-u^t.  The  present  analysis 

does  not  account  for  the  elasticity  of  either  the  cup 
or  the  supports,  largely  due  to  the  lack  of  knowledge 
to  quantify  such. 

In  all  but  one  of  the  present  four  series  of 
calculations,  account  is  taken  of  the  possibility 
of  lack  of  contact  of  the  specimen  with  either  the 
Cup  or  the  supports  (i.e.,  the  tup  and  supports  can 
"push"  the  specimen  but  not  “pull")  at  various  in¬ 
stants  of  time,  as  and  when  the  analysis  may  naturally 
dictate.  In  one  case,  to  study  the  effect  of  the 
above  contact/no-contact  conditions,  the  specimen  was 
held  "fixed"  (i.e.,  the  tup  and  supports  are  always 
in  contact  with  the  specimen). 

Also,  it  is  to  be  understood  that  the  present 
series  of  computations  arc  the  so-called  "simulation" 
or  "generation"  studies  In  the  sense  defined  in  (1|, 
as  opposed  to  tile  "propagation"  studies  performed  in 
( 1 1 .  To  the  uninitiated,  these  terras  imply:  "gener¬ 
ation  study"  means  solving  for  Che  dynamic  stress- 
intensiry  tactor  for  a  crack  which  is  "sude"  to  propa¬ 
gate  with  the  experimenatlly  measured  crack- length 
(or  crack-velocity)  versus  time  history,  while  "propa¬ 
gation  study"  implies  soving  for  the  crack-length 
(or  velocity)  versus  cime  history  using  a  "hypothe¬ 
sized"  dynamic  fracture  toughness  (which  may  or  may 
not  depend  on  crack— salacity  (l))  value. 

The  present  series  of  romput.it ions  are  sum¬ 
marized  in  Table  l. 


Table  I:  "Generation  Studies” 


X.  Data 

Study  1 

No 

Notch-root 

Diameter 

Initiation 

Time 

Boundary 

Conditions 

DTT  1 

0.064mm 

95  user. 

Flxsd 

DTT  2 

0.064mm 

95  usee. 

Contact/ 

No-contact 

DTT  3 

0.064mm 

35  usee. 

Contact/ 

No-contact 

DTT  4 

0.000mm 

92.24  usee. 

Contact/ 

No-contact 

The  crack-length  versus  time  histories  for 
the  above  four  cases  are  shown  in  Fig.  2.  For  the 
cases  DTT  1  and  DTT  2,  the  Aa  vs  t  curves  used  are 
the  same  as  that  in  Fig.  4  of  Ref.  [1],  except  that 
the  boundary  conditions  are  different  as  in  Table  1. 
In  Fig.  4  of  Ref.  [1],  one  date  point  indicating  a 
crack-growth  of  -  2.5mm  at  t-45  usee,  is  Included, 
but  this  point  was  omitted  in  [1]  in  plotting  the 
Aa  vs  t  curve.  In  the  present  study  designated  as 
DTT  3  above,  this  data  point  was  included  in  the 
Aa  vs  t  curve,  and  further,  the  time  of  initiation 
of  propagation  was  chosen  such  that  the  K  value  at 

initiation  -  65MNm  this  curve  is  shown  in  Fig. 

2.  DTT  4  Indicates  the  data  obtained  from  M.F. 
Kannlnen  [4]  for  a  fatigue  pre-cracked  specimen; 
however,  the  initiation  time  was  determined  to  be 
92.24  sec.  by  extrapolation  of  the  supplied  (4] 
experimental  data. 

Prior  to  the  presentation  of  the  results, 
we  Indicate  briefly  the  analysis  procedure.  As  noted 
earlier,  the  present  analysis  of  dynamic  crack  prop- 
atation  is  based  on  the  procedure  developed  by  the 
authors,  and  detailed  elsewhere  [2,  3).  To  supple¬ 
ment  the  mathematical  procedures  In  [2,  31  for  the 
present  case,  we  consider  some  details  of  Imposing 
"contact/no-contact"  boundary  conditions  on  the 
specimen. 

We  designate  the  force  with  which  either  the 
tup  or  the  supports  "push"  the  specimen  as  being 
(+  ve).  Using  the  standard  notation,  the  reaction 
forces  at  the  points  where  displacements  are  pre¬ 
scribed  are  calculated  by: 


P  “  Kq  +  mq .  ( 1 ) 

-  . '  .  ~ 

The  displacement  u  and  reaction  force  P  in  the  time 
step  (n+l)  'are  predicted  by: 
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It  is  noted  that  we  may  use  At  .  ,»4t  *At.  Assume 

n+l  n 

that  P  ,  and  P  are  (+  ve)  and  that  P  ,»P  .  The 
n-l  n  n-l  n 

no-contact  condition  during  the  time  Increment  (n)  • 
(n+l)  is  predicted  to  occur  after  the  sub-increment 

of  time: 
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If  0<At..<At,  during  Che  (n+1)  seep,  we  change  Ac  co 
c 

Ac^  and  perform  Che  analysis  wich  condicion  of  concacC 
and  during  (n+2)  seep,  we  change  Ac  co  Acp  (Atg+Atp 

■Ac),  and  perform  Che  analysis  wich  che  condicion  of 
no-concacc.  This  process  is  repealed. 

An  analogous  scheme  is  used  co  predlcc  Che  tran- 
sicion  from  a  "no-concacc"  co  "concacC"  condicion; 
however,  chis  Clme  by  monicoring  che  dlsplacemencs 
of  che  respeccive  poincs  of  che  specimen  relacive  co 
elcher  che  supports  or  Che  cup. 

RESULTS 


First  we  consider  che  DTT  1  case.  In  Chis  case, 
Che  displacement,  velocity,  and  acceleration  at  the 
point  L  (see  Fig.  1)  are  prescribed.  L^,  and  Lj 

as  shown  in  Fig.  3  are  the  times  when  the  reaction 
force  at  the  cup  becomes  zero  (note  "+  ve  reaction" 
Implies  that  the  tup  is  pushing  the  specimen).  Nega¬ 
tive  "reaction  force"  is  observed  during  times 
Lj<  t<  ,  and  t>L^.  This  phenomenon  can  also  be  ob¬ 
served  in  che  experiments  of  Mall  et  al  [5].  Since 
the  cup  contact  time  measured  in  Che  experlmenC  [1] 
was  about  180  usee.,  the  computation  was  stopped  at 
around  chis  clme.  Fig.  A  shows  che  variation  of  che 
computed  stress  intensity  factors  in  che  present 
"generation"  or  "simulation"  study.  The  times  marked 
by  'v^‘  and  'v^'  in  Fig.  4  are  the  times  when  the 

crack  propagates  with  constant  velocities  v^  and  v^ 

respectively  (see  [1]).  The  Initiation  toughness 
Kjj  (the  terminology  is  summarized  later  in  this 

paper)  obtained  in  this  computation,  as  seen  from 

Fig.  4  is  about  106MNm  ’  .  The  computed  variation 
—  of  input,  scraln,  kinetic,  and  fracture  energies 
wich  clme,  are  shown  in  Fig.  5.  It  is  noted  thac 
In  che  present  procedure  [2,  3]  the  dynamic  K- factors 
are  solved  for,  direedy.  From  this,  the  energy  re¬ 
lease-rate  is  calculated.  Alternatively,  fracture 
energy  is  also  calculated  dlreccly  froR  a  crack-tip 
integral  of  work  done  in  separation  of  crack-faces. 
These  two  procedures,  discussed  in  (6],  were  noted 
to  give  almost  identical  results  for  the  present 
cases).  Ic  is  noted  that  each  of  the  four  energies, 
input,  scraln.  kinetic,  and  fracture,  are  calculated 
Independently,  in  the  DTT  1  case  as  well  as  the  other 
three  cases,  in  the  present  work.  That  these  energies 
"balance"  is  an  independent,  a  posteriori  check  on 
the  present  calculations.  It  is  seen  from  Fig.  5 
that  during  che  periods  L^td^  and  t'Lj  input  energy 

appears  to  actually  decrease,  due  to  "negative  reac¬ 
tion  forces"  during  these  times  as  discussed  earlier. 
We  will  coramenc  further  on  those  energy  variations 
later.  Fig.  6  shows  the  crack  opening  profiles  .it 
various  times.  As  seen,  the  profiles  are  nearly 
linear  except  very  near  the  crack-tip.  This  suggests 
the  possibility  ui  developing  measurement  tech¬ 
niques  from  crack-opening  displacements,  in  dynamic 
tear  testing.  Fig.  7  shows  the  contours  of  equivalent 
stress  oe  in  the  presently  used  singular  element, 

which  is  shown  hatched  in  Fig.  1.  As  seen,  the  stress 
level  reaches  that  of  yield,  a  ,  only  within  the  con¬ 
tour  shown  by  a  broken  line  In  Pig.  7.  The  maximum 
distance  of  this  cotour  is  rough Ly  0.3mm  from  the 


crack-tip;  thus  indicating  the  insignificant  rote  of 
plasticity  in  the  present  problem.  Contours  of 
principal-stress  difference  (which  may  be  used  in  a 
Tresca-type  yield  condition)  shown  in  Fig.  8,  can  be 
observed  to  be  more  or  less  similar  to  those  of 
equivalent  stress  (which  may  be  used  in  a  Mlses-cv^e 
yield  condition)  shown  in  Fig.  7. 

Ue  now  consider  the  case  labeled  DTT  2  in 
Table  1.  It  is  noted  that  the  condition  of  contact/ 
no-contact  was  invoked  in  this  case.  As  seen  from 
Fig.  9,  the  specimen  is  not  in  contact  with  the  tup 
during  the  periods  LjStsL^;  Lj<t<L^,  and  t>Lj  as 

marked.  Also,  it  can  be  seen  from  Fig.  10  that  the 
specimen  is  not  in  contact  with  the  supports  during 
the  times  S^t^S^,  and  t>Sj  a®  marked  in  Fig.  10. 

Comparing  Figs.  9  and  10  ic  is  seen  that  the  maximum 
reaction  force  Pg ,  at  the  support  point,  is  very 

close  to  che  maximum  tup  load,  P^.  The  variation  of 

the  computed  dynamic  K-factor  is  shown  in  Fig.  11. 

It  is  seen  that  the  initiation  toughness,  K_,  is 
-1  5  Id 

again  about  108MNm  ‘  .  However,  prior  to  initiation, 

K  value  appears  to  reach  122MNm  l’^(>K_.)  at  t»82 
1  la 

Usee.  Comparing  Figs.  9  and  10,  it  is  seen  that 
during  the  times  L^xtcS^,  S3<t<L4,  and  t>Lj,  the 

specimen  is  not  in  contact  with  either  che  Cup  or  the 
support,  i.e.,  the  specimen  is  a  free-flying  object! 
Fig.  12  shows  the  variation  of  the  four  eneigy 
quantities;  Input,  kinetic,  elastic,  and  fracture. 
During  LjCtcL^,  Lj<t<L^  and  t»Lj,  since  the  specimen 

looses  contact  with  the  tup,  no  Increase  in  input 
energy  occurs,  as  seen  from  Fig.  12.  It  is  noted 
from  Fig.  12  that  the  total  work  done  at  «180  usee, 
was  about  53  Joules.  This  is  less  than  half  of  the 
experimentally  measured  absorbed-energy  value  of  130 
Joules.  However,  in  the  present  analysis,  no  account 
Is  taken  of  energy  dissipated,  if  any,  in  the  supports 
or  the  tup.  This  discrepancy,  while  of  no  consequence 
in  the  present  "generation-type”  study  (wherein  energy- 
balance  calculations  are  "by-products"  of  the  analy¬ 
sis),  can  have  deleterious  consequences  on  a  "prop- 
atation"-type  study  as  in  [1]  wherein  the  experimen¬ 
tally  measured  absorbed-energy  is  used  in  inferring 
the  fracture  energy,  from  which  a  "plausible"  frac¬ 
ture-toughness  value  is  calculated  [11.  This  is  an 
important  point  to  remember  in  romparing  the  present 
results  and  those  in  [1],  In  Fig.  13,  the  crack¬ 
opening  displacements  in  cases  DTT  1  and  DTT  2  are 
compared.  It  is  seen  that  the  COD  in  DTT  2  oscil¬ 
lates  around  the  values  in  DTT  l  case. 

Now,  we  consider  the  DTT  3  case.  The  tup 
contact-force  variation  is  shown  in  Fig.  14.  Com¬ 
paring  Figs.  9  and  14,  it  is  seen  chat  during 
L^<t<I.7,  the  separation  between  the  specimen  and  the 

tup,  (u^-u^)  b(88er  than  in  the  DTT  2  case.  This 

is  attributed  to  the  higher  compliance  of  the  DTT  3 
specimen  due  to  the  fact  that  growth  initiation  oc¬ 
curs  much  earlier.  Also  the  second  loss  of  contact 
of  the  specimen  and  the  tup  (L,«.t<L.)  and  the  second 

1  a 

peak  of  P  (L  <t«L  )  are  smaller  than  those  In  the 

M  J 

case  of  the  DTT  2  specimen.  Fig.  15  shows  the  dis¬ 
placement  of  the  specimen  from  the  support  point, 
and  the  support  reaction  force.  Comparing  Figs.  10 

and  15  It  is  seen  thac  the  periods  S.-S,  and  S,-S 

1  Z  t.  i 

are  longer  than  those  of  the  DTT  2  specimen:  and  the 

peak  value  of  P  Is  smaller  than  that  of  DTT  2.  Again, 
s 


3 


r 


these  tendencies  can  he  attributed  to  the  earlier 
crack  initiation  in  this  specimen.  Fig.  16  shows  the 
K— factor  variation.  Note  that,  as  shown  in  Table  1, 
the  crack  growth  initiation  time  was  chosen  to  be 

t“35  usee,  such  that  Kr  *65MNm  In  spite  of  this. 

La 

the  K-faetor  variation  in  DTT  3  is  more  or  less  iden¬ 
tical  to  that  in  DTT  2  (Fig.  11)  until  about  t«95usec. 
During  the  period  of  t -95  usee,  to  146  usee.,  the 
K  value  in  DTT  3  decreases  while  chat  of  DTT  2  in¬ 
creases.  The  variation  of  total  work,  strain  energy, 
and  ft?  are  energy  are  shown  in  Fig.  17.  It  is  seen 
that  during  v^<t<v^,  a  very  small  amount  of  energy 

is  consumed  in  the  fracture  process. 

Finally,  we  consider  Che  DTT  4  case.  The  tup 
contact-force  variation  is  shown  in  Fig.  18.  As  seen 
from  Fig.  2,  the  crack  propagates  faster  in  DTT  4 
case  than  in  the  DTT  2  case.  Due  to  this  reason,  the 
changed  compliances,  the  second  peak  of  (L^tcL^) 

is  smaller,  and  the  second  separation  between  the 
specimen  and  the  tup,  l.e.,  (u^-u^  between(L^<t<L^) , 

is  larger  in  the  DTT  4  case  than  in  DTT  2.  The  sup¬ 
port  reactions,  and  separation  between  supports  and 
specimen,  are  shown  in  Fig.  19.  These  are  more  or 
less  similar  to  chose  in  the  DTT  2  specimen.  From 
the  K-faccor  variation  shown  in  Fig.  20,  it  is  seen 

that  the  KTJ  value  Is  about  lllMNm  After 

JLd 

initiation,  K  drops  significantly.  The  energy 
variation  plots  are  given  in  Fig.  21,  from  which  it 
is  seen  thac  Che  total  energy  to  the  specimen  is 
lower  than  in  the  other  three  cases. 

In  the  following  we  attempt  to  draw  some  rational 
conclusions  from  the  above  presented  numerical  data. 

Comments  on  Fracture  in  Impact  Specimens 

At  the  outset,  we  define  the  following  nomen¬ 
clature: 

K^:  Plane-strain  fracture  toughness  under  quasi- 

static  loading 

Kid:  (Variable)  Dynamic  fracture  toughness  for 

a  propagating  crack 

K_d :  Initiation  fracture  toughness  under  dynamic 

loading 

:  "Apparent"  fracture  toughness  at  initiation 

with  a  blunt  notch  under  quasi-static 
loading 

Kqj:  "Apparent"  fracture  toughness  at  initiation 

with  a  blunt  notch  under  dynamic  loading 

KID,min:  Lower  bound-  if  anF-  for  KID- 

A  typical  versus  crnck-ve loc l tv  curve  widely 
reported  in  literature  is  given  in  Fig.  21.  F«u*  most 

brittle  materials,  K  appears  to  be  lower  than 

v  .  ID, min  1  1 


K 


ID, min 


0.5  tu  1.0  K[C. 


(5) 


The  governing  equation  tor  clastodyn.imic  crack  propa¬ 
gation,  currently  considered  co  be  valid  under  quasi- 
static  loading  conditions,  can  be  expressed  bv: 


Kj(t.v)  -  Kid(v) 


(6) 


where  Kjd  is  a  function  of  crack-velocity.  In  general, 

the  above  equation  Is  inconsistent,  from  a  functional 
view-point,  at  the  point  of  crack-growth  initiation. 
For  example,  if  the  "crack"  starts  propagating  (l.e., 
attains  a  finite  velocity  in  zero-time)  from  a  blunt 
notch  with  an  intensity  value  Kq,  we  have: 


K 


1 


Kp(b)  4  K 


ID, min 


where  b  is  the  notch-root  diameter.  Moreover,  for 
initiation  under  quasi-static  loading  from  a  sharp 
notch  or  a  fatigue  crack  [b*0],  wherein  the  crack 
is  assumed  to  attain  a  finite  velocity  in  zero  time, 
we  have: 


Kt  -  K  4  K  , 

I  IC  ID, min 


(8) 


except  in  the  special  when  K1D 

Now,  under  conditions  of  Impact  loading, 

the  initiation  fracture  toughness  K  ,  as  reported 

Id 

in  most  of  the  literature,  is  generally  lower  than 
K^,.  The  values  of  K  as  reported  in  literature 

are  summarized  in  Table  2. 


Table  2:  Initiation  Toughness 
_ Under  Dynamic  Loading 


Investigators 

Material 

Temp  °C 

Loed  Rate 
k  [Mta-1,5/ 

sec.  ] 

Kid/Kic; 

Shabbits  [7] 

A533B 

10 

104 

-0.67 

Kalthoff 
et  al  [8] 

Araldite  B 

Ambient 

104 

'0.73 

i 

Mall  et  al 
[5] 

Poly¬ 

carbonate 

Ambient 

6xl03 

-0.65 

Ireland  [91 

A533 

10 

io6 

-0.63 

Ireland  [9] 

4340  Steel 

— 

20 

Pre-cracked 

Charpy 

-1.00 

Influence  of  the  rate  of  loading  on  is  shown  in 

Fig.  23,  whicli  is  taken  from  Refs.  [7,  9].  As  seen 
from  Fig.  23,  for  most  of  the  materials,  de¬ 
creases  with  increasing  K^.  Fig.  24,  taken  from  Ref. 

[10)  shows  the  comparison  of  K, ,  K  (subscript  I 

d  D 

omitted  due  to  the  fact  that  plane-strain  conditions 
were  not  validly  met),  and  data  for  A533B  steel. 

From  table  2  and  Figs.  23  and  24,  one  may  observe 
that : 


Kld  *  KId(V  ■  °-5  '  U0  KIC-  (9) 


Oner  again,  Kq.  (ft)  1$  functionally  Inconsistent  at 
growth-init lot  ion  associated  with  an  impact  loading. 
Furchrc,  tor  growth-initlaion  from  a  blunt  notch 
und«» r  impact  Loading,  we  have: 


Vb,(VKId  Jnd  V  *  KID.mm«l(n 


However,  utter  a  certain  amount  of  crack  propagation. 
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n)  njv  be  applicable,  In  .1  rate-insensitive  mu- 
r:il.  even  under  impact  loading,  although  Eq .  (6) 

\  expressed  for  quasi-static  loading.  Contrary  to 
'bis,  for  a  rate-sensitive  material  under  Impact 
ie.iJing,  che  governing  equation  for  dynamic  crack 
-r  ipagation  may  differ  from  Eq.  (6),  and  may  be  re- 
-ritcon,  possibly,  as 

Kjft.v)  -  Kid(v,c)  or  KID(v,K)  (11) 

wnere  c  is  the  strain-rate  in  the  vicinity  of  the 
crack-tip  (see  Fig.  25).  If  t  (or  k)  is  small,  or 
the  material  is  almost  rate-insensitive,  the  situations 
in  Fig.  2S  a  and  b  become  practically  the  same. 

It  should  be  noted  that  all  the  analyses  la  Ref. 
(1)  are  of  the  "propagation"  type  in  the  sense  de¬ 
fined  earlier.  In  [1],  the  K  versus  v  curve  for 
4340  steel  is  firsc  assumed  asT 

K1d  «  65  +  0.044  v  (12) 

where  R. „  is  in  MNm  ^ ‘ and  v  is  in  m/sec.  In  [lj, 

the  range  of  crack-velocity  for  which  Eq.  (12)  is 
valid  is  not  indicated.  However,  in  general,  there 
should  be  an  upper-limit  value  for  velocity  as  shown 
in  Fig.  22.  Additional  evidence  of  an  upper  limiting 
value  for  veloclcy  is  shown  in  Fig.  26,  taken  from 

Ref.  (101.  From  Eq.  (12),  KT„  .  -  65MNm“1,5  while 

10, min 

K„  appears  to  be  about  50  from  the  Damage  Tolerant 

Design  Handbook  [  1 1 1 ,  and  Fig.  26.  Considering  the 
near  rate-insensiti  ve  behavior  of  4340  steel,  it  then 

appears  that  K  sK  • K  ,  i50-65MNm 

id  lu  iu,min 

In  order  to  sustain  a  running  cracfc,  the  sur¬ 
rounding  elastic  field  must  produce  plastic  strains 
continually  near  the  advancing  crack-tip  adequate 
for  the  separatlonal  process  (region  R  of  Fig.  25). 

If  K,  .  is  considerably  lower  than  Krn  and  the 
Id  '  ID, min 

elastic  field  near  the  crack-tip  cannot  continually 
supply  energy  to  the  process  zone,  the  crack  will  be 
arrested.  In  Che  above  case  if,  on  the  other  hand, 
che  elastic  field  can  continually  supply  energy  to 
Che  process  zone,  the  crack  will  propagate  with  a 
somewhat  low  velocity  in  the  earlier  stages  of  crack- 
propagation.  This  appears  to  be  verified  in  the 
present  set  of  simulations.  Comparing  che  varia¬ 
tions  in  DTT  2  and  DTT  3  specimens  (Figs.  11  and  16 
respectively)  it  is  seen  that  che  variations  are 

essentially  similar  eventhough  the  crack  in  the  0TT 
3  case  has  been  assumed  to  propagate  with  an  lnltla- 
„l  .5 

tion  value  of  K  *65MNm  .  This  can  be  attributed 
to  the  nacure  of  the  presently  considered  experiment 
in  which  the  elastic  field  near  the  crack  is  always 
growing  in  intensity  during  the  period  t«10  to  80 
usee. 

However,  in  the  DTT  4  specimen,  the  present 
simulation  indicates  Initiation  Choughness  of 

-  HOMNm  (see  Fig.  20).  However,  as  mentioned 
earlier,  the  inici  \t  ion  time  for  che  OFT  4  case  vhlih 
was  unknown  in  che  experiment,  was  determined  by 
extrapolation  of  the  crack-growth  versus  time  data. 
Thus,  as  .in  be  rationalized  from  a  comp. irt. son  of 
the  DTT  2  and  DTT  )  cases,  one  mav  surmise  thar  crick- 
growth  initiation  nav  have  oc.-ur.-d  much  before  t  **»  2 
usee.  (as  used  in  ilic  present  calcul.it  Ions,  see  Fig. 
20).  even  It  with  1  verv  small  velocity.  Dior.  tore. 

■  me  T..1V  postulate  that  even  lor  the  DTT  '•  case.  K 


value  may  actually  have  been  less  than  or  equal  to 

65MNm  This  hypothesis  may  in  fact  be  considered 

to  be  also  supported  by  the  results  of  [1].  Refferlng 
to  Fig.  8  of  (1),  it  is  seen  tha  in  the  "propagation" 
calculation,  the  use  of  the  relation  Kjp*65+0 ,04(^v 

actually  gives  a  much  better  agreement  with  the  experi¬ 
ment  in  terms  of  crack-initiation  time  rather  than  the 

use  of  Kjp*170KNm  Note  that  the  experimental 

data  for  the  first  measurement  of  crack  grwoth  is 
missing  in  Fig.  8  in  Ref.  (1).  Contrary  to  this, 
after  a  certain  amount  of  crack-growth,  use  of 

Kjjj*170MNm  appears  to  give  a  good  agreement  of 

the  "propagation"  calculation  (1)  with  the  experimen¬ 
tal  data. 

However,  as  noted  earlier,  the  value  of  Kt 
-1  5 

•170MNm  ‘  was  "derived"  in  (1)  from  an  energy  bal¬ 

ance  consideration:  the  total  energy  absorbed  during 
the  impact  fracture  of  specimen  was  130  Joules;  the 
energy  imparted  to  an  "intact”  unconstrained  specimen 
from  an  elastic  collision  between  the  top  and  speci¬ 
men  was  79  Joules;  and  the  difference  (130-79)-51 
Joules  was  assumed  to  be  consumed  as  fracture  energy. 
From  this  assumed  fracture  energy  and  average  frac¬ 
ture  toughness  value  of  l70MNm  was  derived  in 
[1].  Thus,  the  propagation  analysis  (1]  represents 
an  a  priori  "energy-balance”  condition,  in  which  other 
sources  of  energy  dissipation  are  Ignored.  On  the 
other  hand,  as  already  noted,  in  the  present  "genera¬ 
tion"  calculations,  an  energy  balance  relation  is 
an  a  posteriori  "by-product"  of  the  calculation  it¬ 
self  in  which  the  experimental  crack-growth  history 
Is  simulated.  In  connection  with  Fig.  12,  (DTT  2 
specimen,  which  is  also  the  basis  for  the  computation 
K1d«170  in  Ref.  (1])  it  is  noted  that  the  computed 

total  input  to  the  specimen  is  53  Joules,  which  is 
less  than  half  of  the  experimentally  measured  value 
of  130  Joules.’  This  brings  to  question  che  negleccion 
of  ocher  dissipated  energies,  and  hence  the  value 

of  KI0-l70MNm  hypothesized  in  (!]. 

The  present  analysis,  and  the  above  dis¬ 
cussion,  appear  to  lead  us  to  believe  that  the  use  of 

K._  versus  v  curve  with  K-,,  ,  *  65MNm  and  K_.< 

ID  .  ID, min  Id- 

65MNm  *  ,  with  che  curve  being  such  that  the  upper 
value  of  v  is  limited  (as  in  Fig.  22),  may  be  war¬ 
ranted.  Suppose  this  "saturated"  curve  is  used  in 
the  "propagation"  calculation;  then  the  crack  will 
start  to  propagate  at  t  *  25  1  35  usee.  During 
most  of  the  period  from  initiation  to  about  t  -  100 
usee.  ,  the  crack  will  propagate  with  .1  relatively  slow 
velocity,  as  discussed  earlier  (l.e.,  crack  grouch 
will  be  small  until  about  t  -  100  usee.).  The  crack 
will  attain  a  highcr  velocity  because  of  the  nature 
of  the  testing  system  considered,  in  chat  is  In¬ 
creased  continually  by  the  dropping  tup.  Then  the 
crack  speed  will  saturate.  in  addition  to  the  prop¬ 
erty  of  the  "saturated"  K^vs  v  curve,  there  may  be 

other  faccors  to  limit  the  crack  speed  in  the  present 
type  of  specimen.  One  of  these  is  the  presence  of  the 
compressive  stress  field  in  the  ligament  of  t lie  speci¬ 
men,  is  can  be  seen  from  the  momentum  balance  condi¬ 
tion. 

Cootrarv  to  the  impact  loading,  mJer  quisi- 
static  loading  tile  analysis  using  ♦  J.  v 

gave  guild  igreement  with  the  experiment  in  (LI.  tn 
this  .  ise,  since  the  tup  was  fixed  to  correspond  to 
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Kg-lOSMNnT1'5,  the  value  of  Kj  during  crack-propa¬ 
gation  under  quasi-static  loading  condition  is  always 
lower  than  108MNm  [Kj(t, v)<K^] .  This  suggests 
that  Kj^-65  *■  0.044  v  is  valid  for  lower  values  of 


velocity,  while  this  linear  relation  may  be  Invalid 
for  larger  velocities,  and  eventually  the  velocity 
will  be  limited,  leading  to  a  satuaraced  versus 
v  curve. 
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Numerical  simulations  of  crack-propagation  histories  in  four  cases  of  dynamic 
tear  test  experiments  on  4340  steel  are  performed.  The  influence  of  the  loss 
of  contact  of  the  specimen  at  various  times  with  either  the  supports  or  the 
tup  or  both  is  critically  examined.  In  each  case,  the  variation  of  the  dynamic 
K-factor,  for  the  simulated  crack-propagation  history,  is  directly  computed. 

The  results  are  discussed  in  the  light  of  current  controversies  surrounding  the 
dynamic  fracture  toughness  properties  governing  crack-propagation  under  impact 
loading.  It  is  concluded  that  these  controversies  may  not  be  fully  warranted. 
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